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We present a microscopic quantum theory of intersubband polarons, quasiparticles originated 
from the coupling between intersubband transitions and longitudinal optical phonons. To this aim 
we develop a second quantized theory taking into account both the Frohlich interaction between 
phonons and intersubband transitions and the Coulomb interaction between the intersubband tran- 
sitions themselves. Our results show that the coupling between the phonons and the intersubband 
transitions is extremely intense, thanks both to the collective nature of the intersubband excitations 
and to the natural tight confinement of optical phonons. Not only the coupling is strong enough 
to spectroscopically resolve the resonant splitting between the modes (strong coupling regime) , but 
it can become comparable to the bare frequency of the excitations (ultrastrong coupling regime). 
We thus predict the possibility to exploit intersubband polarons both for applied optoelectronic re- 
search, where a precise control of the phonon resonances is needed, and also to observe fundamental 
quantum vacuum physics, typical of the ultrastrong coupling regime. 



I. INTRODUCTION 

The theory of polarons, the quasiparticles describing 
electrons in a polarizable medium, dates back to the early 
days of quantum theory^, and it has been an active field 
of research ever since^. 

In this paper we will develop a microscopic theory of 
intersubband polarons, that is, a theory of intersubband 
transitions coupled to longitudinal optical (LO) phonons 
in semiconductor quantum wells. 

The coupling between intersubband transitions and 
LO-phonons is relevant for a number of optoelectronic ap- 
plications, as it determines the lifetime of carriers in ex- 
cited subbands^. In particular a precise knowledge of LO- 
phonons intersubband scattering rates is important in 
the engineering of heterostructures for quantum cascade 
lasers'*. Normally optoelectronic devices are designed to 
avoid being in resonance with optical phonon transitions, 
due to the high absorption between transverse and longi- 
tudinal optical phonon frequencies (Restrahlen band) . A 
notable exception is provided by quantum cascade lasers 
operating near such optical resonances^, in which in- 
stead the transitions between different subbands are al- 
most resonant with LO-phonon modes. 

Even if the coupling between intersubband transitions 
and LO-phonons in semiconductor quantum wells has in- 
deed received some attention^^^ and intersubband po- 
laron resonances have been clearly and unambiguously 
observectiS, to the best of our knowledge, there is no mi- 
croscopic theory of such excitations, as the spectra of 
intersubband polarons are normally calculated with in- 
direct methods. While such methods allow to calculate, 
at least qualitatively, the polaron dispersions, missing a 
microscopic description makes it difficult to study more 
complex phenomena as nonequilibrium physics, quantum 
vacuum effects or quantum phase transitions. 

Using a second quantization formalism, we will re- 
duce the full electron-phonon Hamiltonian to a quadratic, 
bosonic form, from which we will then calculate the po- 
laron dispersions. In order to accomplish this task, we 



show that Coulomb interaction between electrons in con- 
duction subbands naturally separates into a dominant 
and a perturbative part, accordingly to the number of 
electrons that can participate to each transition. 

Moreover we will show how, thanks to the tight con- 
finement of LO-phonons, the coupling between intersub- 
band transitions and phonons can easily be in the ultra- 
strong coupling regime, a regime characterized by a cou- 
pling strength comparable to the bare frequency of the 
excitations^. Such fact can have interesting observable 
consequences, as a whole new range of physics is a pri- 
ori observable in this regime: spectral deformations^^^— , 
quantum vacuum emission phenomenai^^— , electrolu- 
minescence enhancemen t ^^d^ and even quantum phase 
transitions--^--. 

This article is organized as follows: in Sec. |TT] we 
will develop the general theory of the coupling between 
intersubband transitions and LO-phonons, that we will 
then apply in Sec. IIIII to the case of GaAs quantum 
wells, showing how the ultrastrong coupling regime can 
be reached even with such relatively weakly polar ma- 
terial. In Sec. IIVI we will compare the calculated dis- 
persions with the ones obtained using an homogeneous 
dielectric function approach. Finally a few considera- 
tions on the impact of our results and on possible future 
developments will be drawn in Sec. [V] 

II. THEORETICAL FRAMEWORK 
A. Superradiant excitations 

In 1954 Dicke^ noticed that a set of coherently excited 
identical dipoles relaxes radiatively much faster than a 
single, isolated one. This is due to the phenomenon of 
superradiance: N identical dipoles behave as a single col- 
lective dipole ^/N times bigger. 

The concept of superradiance has been thoroughly ap- 
plied to the study of intersubband polaritons^ir— in mi- 
crocavity embedded quantum wells. In such systems the 



light couples to a collective electronic excitation and, as 
expectedii, the strength of the coupling between light 
and matter is proportional to the square root of the num- 
ber of electrons involved. 

We will study the coupling of intersubband transitions 
with longitudinal optical phonons, considering also the 
role of Coulomb electron-electron interaction. Such cou- 
plings are extremely rich and, in order to limit the com- 
plexity of our investigation, we will need to determine 
which scattering channels are dominant and which are 
negligible. In general, if N electrons undergo a certain 
transition in a coherent way, the strength of the cou- 
pling is enhanced by a factor -x/ZV- Transitions involving 
a macroscopic number of electrons will thus be strongly 
enhanced and, for this reason, they will be treated ex- 
actly within an Hamiltonian formalism, while the others 
(involving only few electrons ) will be treated perturba- 
tively (or ignored altogether). The degree of collective 
enhancement of a scattering process will be evaluated 
looking at the number of electrons that can undergo the 
transition given fixed amounts of transferred impulsion 
and energy. In Fig. [T]we show a few illustrative exam- 
ples of collective and non-collective transitions. 

In the case of the Coulomb interaction, as we will see, 
only the intersubband terms, responsible for the depolar- 
ization shift, are collective. We will thus treat them ex- 
actly in the Hamiltonian, while considering all the other 
terms in a RPA linear response approach. 



B. Free fields 

We will consider a symmetric quantum well of length 
Lqw in a bulk of height Lbk, S will be the surface of the 
sample. For the moment, we will limit ourselves to the 
case of a single quantum well, the general case of multiple 
wells will be addressed later in this Section. 

The quantum well is supposed to be doped in such 
a way that its Fermi level is between the first and the 
second conduction subbands, separated between them by 
the intersubband gap energy huji2. 

We will develop our theory using a zero temperature 
formalism (T = 0), anyway our results will remain quan- 
titatively accurate while the thermal population of the 
second subband remains negligeable. Depending on ma- 
terial parameters and doping level, this could imply the 
necessity to perform experiments in different kinds of 
cryogenic environments. 

Electron states will be indexed by the subband index 
J and by the value of the in-plane wavevector k. Their 
wavefunctions will be given by 

tpj,k{p,z)^Xj{z)—^, J = 1,2, (1) 

where, for simplicity, we will choose Xji^) to be real and, 
due to the symmetry of the quantum well, the Xj (z) have 
well defined and opposite symmetry. Wavefunctions in 
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FIG. 1: a) An example of collective transition; an intersub- 
band transition with small transferred momentum. All the 
electrons can undergo a transition resonant approximately at 
the same energy, b) An example of non-collective transition: 
an intrasubband transition with transferred momentum much 
smaller than the fermi wavevector hkp. The majority of elec- 
trons, being Pauli blocked, do not participate to the process, 
c) Another example of non-collective transition: an intrasub- 
band transition with transferred momentum of the order of 
2hkF, that is the minimum to allow all the electrons to un- 
dergo the transition avoiding Pauli blocking. Anyway there 
is a large energy spread between the different single electron 
transitions. For an infinite potential well, electrons on the two 
opposite borders of the Fermi sea, with initial momenta zthkp 
parallel to the transferred momentum, have initially the same 

energy = Ef, where m is the electron effective mass 

and Ef the Fermi energy. After the transition they will end 
up with final momenta hkp and 3hkF, corresponding to final 

^2^2 9^2^2 

energies — Ef and = 9Ef respectively. This im- 

plies that, even if the transition is not blocked, only a small 
fraction of the single electron transitions can be resonant at 
the same time. In both cases, given that only few electrons 
can participate to the collective transition, the superradiant 
enhancement factor will be small. 



Eq. ([TJ are chosen as basis for second quantization, the 
creation operator for an electron in the state described 
by Eq. ^ will be denoted as c^y^. The free Hamiltonian 
of the electron gas in the two considered subbands thus 
reads 

Hei = ^'^j(k)4kCj,k, (2) 

where wi(k) — and W2(k) = a;i(k) -I- wi2. In Eq. 
(12), as well as in the rest of this paper, we will omit the 
electron spin index. This is justified by the fact that 
all interactions we consider are spin conserving. Given 
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that we will consider only in-plane wavevector exchanges 
q much smaller than the typical electron wavevector k, 
we can make the approximation u)j{\i + q) ~ ujj(k.)^ and 
introduce the operators describing intersubband transi- 
tions with a well defined and dispersionless energy Hlu 
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(3) 



C. Electron phonon interaction 

Interaction between electrons and LO-phonons can be 
described using the Frohlich Hamiltonian^ 



where 



-i(qp+g2z) 



2eoepSLBK ^ y/q'^ + q 



-dl^^+h.c, (8) 



(9) 



where N is the number of electrons in the quantum well. 
In the dilute regime, that is if the number of excitations 
is much smaller than N, the b]^ operators are bosonic 



3q,q'- 



(4) 



At higher excitation densities, that are out of the scope of 
the present work, saturation effects start to appear and 
corrections to Eq. Q have to be taken into account. For 
a detailed analysis of nonbosonicity effects in intersub- 
band transitions we invite the interested readers to refer 
to Ref. 
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Using Eq. ([3]) we can rewrite the Hamiltonian of the 
free electron gas in Eq. ^ in terms of bosonic intersub- 
band excitations 



(5) 



In this work we are interested in the resonant case 
in which a;i2 is equal, or close, to the LO-phonon fre- 
quency ujlo- We can thus neglect confinement ef- 
fects on the phonons and consider bulk values for their 
frequencie a^^'^" . We will thus describe LO-phonons by 
means of the three dimensional boson operators dq^q^ 



(6) 



indexed by their in-plane and out-of-plane wavevectors. 
While we know that LO-phonon modes are confined in- 
side the quantum well, we do not need to impose this 
constraint in the mode definition because, as we will see, 
intersubband transitions end up coupling with linear su- 
perpositions of phonon modes that are anyway confined 
inside the quantum well. We will consider only the case 
of one single longitudinal optical branch, the expansion 
to the case of multiple branches not presenting any fun- 
damental difficulty. 

Moreover, we are interested only in phonons with small 
in-plane wavevectors (in order to couple with coherent in- 
tersubband excitations), we can thus ignore phonon dis- 
persion and write the free phonon Hamiltonian as 



ph 



nUJL0d\^q^dq^q^ 



(7) 



where and are respectively the static and high fre- 
quency dielectric constants^. 

The Hamiltonian in Eq. ([8]) can be written in sec- 
ond quantization (neglecting incoherent intrasubband 
scattering^) as 



E 



F{q. 



(C^q,g, + rf-q -?J(cl^kC2,k-Hq + 4,k-qCl.k), 

where we have defined 



(10) 



dzxi{z)x2{z)e~^'' 



(11) 



From Eq. pO|) we see that, due to the three dimen- 
sional character of the LO-phonons^°, each electronic 
transition couples to multiple phonon modes, indexed by 
different values of the wavevector along the growth direc- 
tion. It is thus convenient to introduce second quantized 
operators corresponding to the particular linear super- 
positions of phonon modes that are coupled to electronic 
transitions 



-y 



Fiq.)dl 



Fiqz)dq 



whose spatial wavefunctions along the z axis are 



-y 



fqiz) 



1 



BK 



E 



F{q.)e 



(12) 



(13) 



From Eqs. and (ITU) we see that the intersubband 

transitions naturally couple to phonon modes localized 
inside the quantum well (it is easy to verify that •^q{z) 
vanishes to the first order in g if z is outside the common 
support of xi and x^)- 

The normalization factor A can be fixed imposing 
bosonic commutation relations for the operators 



-y 



F{q.)F{-q'^] 



A 



5q,q'LBKl{q) 



2Aq 



(14) 



and thus 



A 



Lbk I{q) 



a) 



where we have defined 

I{q)= I dzdz'xi{z)X2{z)X2{z')xi{^')e-'^'-'-'^- 



(15) 



(16) 



We can thus write Hamiltonians in Eq. ([7]) and Eq. 
(1101) in terms of coherent and b^^ operators as 



b) 



(17) 



and 

q 

where 



2DEGfn^LO 



4ene 



N 

J^2DEG — -g-i 



(&]j + 6_q)(rq+rlq), 



(18) 



is the density of the two dimensional electron gas. In 
order to pass from Eq. (O to Eq. (IT7|) . we are neglecting 
higher order phonon modes confined inside the quantum 
well. This is justified by the fact that we limit ourselves 
to long-wavelength modes. 



D. Superradiant Electron-Electron interaction 



i, k + q 



irimnj 




FIG. 2: a) Index convention of the matrix element Vq"^"-' . 
Two electrons in subbands j and n, with momenta k and k' 
are scattered into subbands i and m, with momenta k + q 
and k' — q respectively, b) Graphical representation of the 
four qualitatively different kinds of scattering processes from 
Eq. (p2l). 



In order to treat the Coulomb electron-electron in- 
teraction we start by the second quantized form of the 
Hamiltonian describing the Coulomb interaction^ (see 
Fig. [5] (a) for a graphical representation of the interac- 
tion coefficients) 

= J E E ^r"'4k+q4,k'-cCn,k'c„k, (19) 



2, J, m,n— 1,2 q,k,k' 



where 



The other elements can be evaluated as 

T/1122 T/1212 T^2121 _ T^2211 

q 1 1 1 



(22) 



ta1221 _ T^2112 
— ''q 



yiiii = — 

2eoeoog 



2eoeoog 



2eoeoog' 
dzdz'xl{z)xl{z')e-'i\^-^'\ 



yirnnj _ 

' 2eoeoog 



dzdz'xi{z)xj{z)Xm{z')Xn{z')e '^''^ \ 

(20) 



2eoeoo'7 



dzdz'xl{z)xl{z')e-^\^-''\, 
dzdz'xl{z)xl{z')e-''\^-^'\ 



is the two dimensional Coulomb matrix element. It is 
important to notice that in Eq. (|20|) we used the high 
frequency dielectric constant too instead of the static one. 
This is due to the fact that es includes the effect of the 
coupling to LO-phonons, that are already treated exactly 
in the Hamiltonian. 

Due to the symmetry of the wavefunctions, a certain 
number of matrix elements in Eq. (j20p can be seen to be 
zero, in particular all the matrix elements with an odd 
number of 1 and 2 indices 



1112 



2111 



1121 



2212 



1211 



2122 



2111 



1222 



0, 
0. 



(21) 



where /(?), defined in Eq. (jl6p . is the same integral we 
encountered studying the electron-phonon Frohlich inter- 
action. 

The four distinct nonzero possible values of the matrix 
elements correspond to different kinds of scattering pro- 
cesses. In Fig. [5] (b) a graphical representation for each 
of these processes is shown. 

It is important at this point to notice a major differ- 
ence between the elements in the first line of Eq. ([^^ 
and the others. The elements in the first line (upper 
left subpanel in Fig. [5] (b)) represent intersubband ex- 
citations: each electron is scattered from one subband 
to the other. Such processes, responsible for the depo- 
larization shift^, describe a superradiant process in the 
sense defined above, that is, at least for small values of 
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q, a great number of electrons can coherently undergo 
the same transition, approximately at the same energy. 
This is not the case for the interactions described in the 
other lines of Eq. (1^^ . that instead describe intrasub- 
band excitations that, either due to Pauli blocking or to 
the non-flat energy dispersions, involve only few electrons 
(see Fig. [T] for a graphical visualization of this crucial 
point). 

Our previous discussion on superradiant processes thus 
implies that the terms in the first line of Eq. (1^^ strongly 



dominate over the others due to their superradiant en- 
hancement. For this reason we have to treat them exactly 
in an Hamiltonian formalism, while we can limit ourselves 
to treat the others within a perturbative approach. 

Here we will thus construct an exact, Hamiltonian ap- 
proach, to treat the effect of the depolarization shift 
terms, neglecting the others. We will analyze later the 
effect of the intrasubband terms. 

Let us start to rewrite the depolarization shift part of 
Eq. ([T^ in a more useful form 



= X! JTT^ (4.k+qcI,k'-qC2.k'C2,k + 4,k+q4,k'-qCl,k'C2.k + 4,k+q4,k'-qC2,k'Cl,k + 4,k+q4,k'-qCl,k'Cl,k) 



q,k,k 



, 4eoeoog 
e^I{q) 



A, , „ (cl,k-f-qC2,ke| i,,_qC2,k' + cl.k+qC2,k4.k'_qCl,k' + C2,k+qCl,kCl,k'-qC2,k' + 4,k+qCl,kC2.k'-qCl,k' 



q,k,k 



E 



Ne^I{q) 



(23) 



We see from Eq. ([231) that, thanks to its collective, super- 
radiant nature, the depolarization shift can be naturally 
written in terms of the bosonic intersubband excitations 
defined in Eq. ([3]) as 



He — ^ N2DEG 



4ene 



(24) 



where we have neglected the last constant term, that sim- 
ply shifts the ground state energy. 



E. Residual Electron-Electron interaction 

We treated exactly the depolarization shift terms of 
the Coulomb interaction in a bosonic excitation formal- 
ism. Moreover, we showed how the terms other that the 
ones responsible for the depolarization shift are strongly 
suppressed, due to their lack of collective enhancement 
and can thus be treated perturbatively. 

Here we will study the perturbative effect of such 
residual Coulomb contributions, due to the intrasubband 
terms in the last three lines of Eq. ((22)) (schematized in 
the last three panels of Fig. [11(b)). 

An important result due to Lee and Galbraith2'^''^'', is 
that such intrasubband terms do not contribute to the 
screening of the intersubband ones at the level of the 
random phase approximation (RPA). This can be seen 
writing the Dyson equation for the dynamically screened 
Coulomb potentiaP^ '^gi'^) 



V™"J'(w) 



^y;"m^'^(w)V™(a;), (25) 



where n*'"(a;) is the RPA polarization function. In the 
case of an intersubband contribution (e.g., V^^^^), Eqs. 
([2T|) and ([22|) imply that 



vl'^^c) = i/"22 + ^y/^^2n-(c.)vf2'^(^) (26) 

rs 

= v,"^^ + ^T^/'-^^n-(c.)vf2'-(a.) 

= yn22 ^ yn22(ni2(^) + nf {u.))V'^''' {u^) . 



We have thus 



V,"^^(c.) 



T/1122 



i-y,"22(ni2(^) + n2i(..))' 



(27) 



from which we see that the intrasubband Coulomb terms 
(V,"",y,2222^ y^i22i ^^^^ y^2ii2) intervene in the 

renormalization of the intersubband terms. 

An analogous reasoning can be done for the phonon- 
electron interaction. Calling Mq^q^ and Mq^q^{uj) the 
bare and screened version of the potential defined in Eq. 
([TU[) . we have the Dyson equation 

Mq,q^ (OJ) = Mq,q^ + ^ y/'-^H^ (..)A^ {uj) (28) 

rs 

= Mq,q^+Vq''^\nf{0j)+nf{u;))Mq^q^{0j), 

and thus the formula for the screened potential is 

■^^''^'^'^^ = i-Vq^^^Hn/{i^) + nf{co))- ^^^^ 

Being the RPA screening only due to terms that are ex- 
actly treated in the Hamiltonian, we can thus neglect the 
screening due to the two dimensional electron gas. 
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F. Hopfield-Bogoliubov Hamiltonian 

Putting together Eqs. ©, (HH) and ([24]) we arrive to 
the full Hamiltonian for the intersubband transitions-LO- 
phonons system 



electronic part of the mixed excitations is an intersub- 
band transition. Hamiltonian in Eq. p2p can thus be 
put in the diagonal form 



i=±.q 



(37) 



(30) 



' N2DEG^LO 



e2 /(g), 



where the pj^q are the annihilation operators for the two 
polaronic branches, given by a linear superposition of 6^, 



[h^ +b ){r^ + r ) ''q' ^-q ^"^"^ '"-q operators and E/^ is the energy of the 



4eoep 9 



+N2DEG 



4eoeoo 9 



The Hamiltonian in Eq. pop can be rewritten in a 
more compact form by introducing the intersubband 
transitions-LO-phonons coupling coefficient and the 
Coulomb coefficient D 



^ — \ N2DEG^LO 



Aeoeph q 



(31) 



D ^ N- 



2DEG 



I{q) 
4eoeoo q 



where we have dropped the dependences over the 
wavevector as we are interested in the long wavelength 
limit (from Eq. we can verify that limg^.o tends 
to a constant value). 

Using Eq. ([nj), Eq. ([301) can be written as 

H = ;i^wi26lj6q + WLO<rq + r!(6]j + 6_q)(rlq + rq) 
q 

+D{bl + b_^){bl^ + b^), (32) 
that can be cast in matrix form as 



(33) 



where the column vector of operators Vq is defined as 

«q= [5q,rq,fe^_q,rlj^, (34) 

rj is the diagonal metric 

?7 = diag[l, 1,-1,-1], (35) 
and the Hopfield-Bogoliubov2S matrix Hq is defined as 



Hq — 





2D 


n 









-n 


-CJ12 - 2D 


-n 





-n 





(36) 



Diagonalizing the matrix in Eq. (1361) will yield the fre- 
quencies of the normal modes of the system uj± , that are 
usually called polarons^i^^. In our case we will name 
them more properly intersubband polarons, because the 



new ground state relative to the one of the uncoupled 
system. 



G. Coupled Ground State 

The coupling between the intersubband transitions and 
the LO-phonons does not modify only the system's reso- 
nances but it also qualitatively modifies the nature of its 
ground state. It is easy to verify that, if |0) is the ground 
state for the uncoupled phonons and intersubband exci- 
tations, defined in the usual way as 



then 



MO) = rq|0)=0. 



Pj,q|0) ^ 0, 



that is, |0) is not the ground state for the coupled system. 
The real ground state of the Hopfield matrix in Eq. (|36p , 
that has been thoroughly studied in Ref. [lH, has the 
form of a two modes squeezed vacuum. 

Still, thanks to the bosonicity of the system, such new 
ground state does not influence the response of the sys- 
tem, that can be described as a gas of free bosonic exci- 
tations (from Eq. ([37| ). A notable exception is the case 
in which the parameters of the system are nonadiabati- 
cally modulated in time. In this case the sudden change 
in the ground stated can have observable effects, like the 
emission of quantum vacuum radiatior^i^ii^. 

It is also interesting to notice that, from Eq. ([5T|). we 
can write the Coulomb coefficient D as 



D = 



> 



(38) 



As it has recently been shown in Ref. Eq. (|55|) im- 
plies that the ground state of the system will not undergo 
a Dicke phase transition, regardless of the strength of the 
coupling. 



H. Multiple Quantum Wells 

Until now we considered the case of a single quantum 
well. This choice has been motivated by the fact that, 
as we will show, the presence of multiple wells does not 
modify our results. 
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Given that we are considering rather large quantum 
wehs (in order for the transition to be resonant with the 
LO-phonon mode), the optical phonon spectrum is not 
modified^*' and the optical phonon modes we consider are 
confined in each quantum well. 

This is a rather important difference between the in- 
tersubband polaron case we consider in this paper and 
the physics of intersubband polaritons. For intersubband 
polaritons, the electromagnetic mode coupled to the in- 
tersubband transitions extends over all the structure. It 
thus couples to all the electrons, regardless of the quan- 
tum well they are in. This means that the only mean- 
ingful parameter for intersubband polaritons is the total 
density of electrons, and the light-matter coupling thus 
scales as nQwN2DEG, where tlqw is the number of 
quantum wells inside the microcavity. 

In the present case instead, being the phonon modes 
confined inside each quantum well, electrons in different 
wells are completely decoupled. This can also be inferred 
from the coupling integral in Eq. (f20| . This integral does 
vanish, at least in the long wavelength limit (first order in 
g), if the wavefunctions for the two integration variables 
z and z' do not have a common support, i.e., if the two 
interacting electrons are in different quantum wells. 

This means that, contrary to the intersubband polari- 
ton case, the intersubband polaron interaction scales only 
as V N2D EG and growing multiple quantum wells in the 
same sample will not increase the coupling. 



III. RESULTS 

In order to obtain some numerical predictions from 
Hamiltonian in Eq. (|32p . we need to fix a few parameters 
concerning the material and the quantum well. 

For sake of simplicity we will consider the quantum well 
to be correctly approximated by a rectangular, infinite 
potential well of length Lqw- We thus have 



hhJi2 = 



(39) 



and the electronic and phononic modes profiles are given 

by 




2 TTZ 

sm( ), 

Lqw Lqw 



(40) 



2 . , 2ttz ^ 
sm(- ), 



16 . TTZ , 

— sm3(- ), 



inside the well and zero outside. As explained in Sec. |TI1 
we see here explicitly that the intersubband transitions 
couple to a linear superposition of phonon modes that 
is localized inside the quantum well (the cubic sinus in 



the third line of Eq. pO)) comes from the integral of the 
first two, as can be verified performing the integral in Eq. 

(HI). 



Inserting Eq. PU)) into Eq. (IT51) and performing the 
integral we have 



lini/(g) ^ ^qLqw 



(41) 



In Fig. [3] we plot the normalized coupling ^j^— as a 
function of the density of the two dimensional electron 
gas, for a GaAs quantum well. In the inset of Fig. [3] 
we instead present a comparison of the values of 
at room temperature, for different semiconductors of tfie 
III-V and II- VI groups'**^, as a function of the respective 
LO-phonon energies, for a reference doping N2DEG = 
IQi^cm-^. 

In Fig. |4]there is a plot of the intersubband polaron fre- 
quencies ui± as a function of the intersubband frequency 
UJ12, in GaAs, for N2DEG = 10"^^cm~^. Notice that, due 
to the effect of Coulomb interaction, the resonant anti- 
crossing is not at LO12 — lulo but at a lower frequency. 
In the inset of the same figure we plot the same quantity 
as a function of the electron density. The length Lqw 
has been chosen in this case to have the two uncoupled 
modes at resonance (wi2 = ljloj that is Lqw — 23nm). 

It is clear from the figures that intersubband polarons 
are not only strongly coupled, having coupling constants 
much larger than their linewidth (usual linewidths being 
not bigger than a few meV), but they are indeed in the 
ultrastrong coupling regime, with values of the normal- 
ized coupling comparable or larger than the best 
ones reported in the literature. For physically realizable 
levels of doping, coupling values of a few tenths of the 
bare frequency of the excitation lulo are predicted in 
GaAs, and it seems that values much larger can be ob- 
tained using more polar materials. The reason of such 
large coupling can be found in the superradiant nature 
of intersubband excitations and in the natural confine- 
ment of the phonons inside the quantum well, that gives 
an extremely small mode volume, when compared with 
what can be obtained with photonic microcavities. 

The consequences of our results can be multiple, both 
for fundamental and applied research. On the fundamen- 
tal side, intersubband polarons could become a new labo- 
ratory to test quantum vacuum physics, typical of the ul- 
trastrong coupling regimei^. On the applied side our the- 
ory can be naturally exploited in the study of quantum 
cascade lasers working in or near the Restrahlen band. 
It can, for example, help explaining the anticrossing ob- 
served in Ref. |5|, near the LO-phonon frequency. More- 
over the capability to strongly modify the LO-phonon 
spectrum could have an impact on the performances of 
optoelectronic devices, as the electron-LO-phonon scat- 
tering rate determines the lifetime of carriers in excited 
subbands^. 
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IV. COMPARISON WITH DIELECTRIC 
FUNCTION THEORY 



In the previous Sections we have developed a detailed 
microscopic theory for the intersubband transitions cou- 
pled to LO-phonons. Here we will compare the disper- 
sions obtained from the microscopic theory with the ones 
obtained with an homogeneous dielectric function theory, 
as the one used in Ref. fl0'|. 

The propagation of an electromagnetic wave in a dis- 



persive, homogeneous medium obeys the equation 

div[i:>(a;)] = dw[e{uj)E{u)] = 0. (42) 

This implies that it is possible to have propagating longi- 
tudinal waves, like polarons, only at frequencies for which 



(43) 



where indicates the real part. 

The z component of the dielectric function of an ho- 
mogeneous medium filled with quantum wells is given 

byia^ 



e(c^) 



where 



^00 9 



42 



LOr 



00-'-'' 



(44) 



(45) 



is the plasma frequency of the two dimensional electron 
gas and di2 is the intersubband dipole 



(ii2 = e / dzxi{z)zx2{z). 



The equation 



thus reads 



3?[e(w)] = 0, 



(46) 



(47) 



(48) 



As Eq. p4|) neglects both the dielectric response in the 
X — y plane and the non-homogeneity in the z direction, 
in order to recover the same result from our microscopic 
approach, we will have to consider only phonon modes 
with q = and -> 0. From Eqs. ([TT]) and (|46]) we 
thus have 



X12 



(49) 



Following exactly the same procedure of Sec. |TI1 but 
with the F(g) defined in Eq. p9|) and considering only 
the 92 mode, we get 




D 



N2DEGdl2 



(50) 



(51) 



In order to obtain the polaronic eigenfrequencies we 
have to diagonalize the matrix in Eq. (I36|) using the 
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coupling coefBcients for the homogeneous hmit defined in 
Eqs. ([50]) and ([5T|) . We thus obtain the secular equation 

Lu^ -lj^{ujIo+uj12 + '^Dlji2) (52) 

+w|oWi2 + 4:DuJi2UjIq - 4ri^lLJi2WLO = 0, 

that, using Eqs. (|5ip and (|45l) can be put into the form 

4 2/2 I 2, 2\i 2 2, 22 ^oo „ 

(53) 

Equating the coefficients of Eqs. (1151) and ([551) . we obtain 
^'to = ^lo'-^, (54) 

that is the well known Lyddane-Sachs- Teller relation^. 
We have thus proved that the homogeneous version of 
our theory gives the same results as the homogeneous 
dielectric function approach. 

It is anyway important to notice that the homogeneous 
limit in not exact, as a quantum well is, by definition, 
spatially inhomogeneous. Ignoring the higher qz modes 
leads to underestimate the intersubband dipole of a factor 
roughly equal to v^. 



V. CONCLUSIONS 

In this paper we have developed a microscopic the- 
ory of intersubband polarons, mixed excitations result- 
ing from the coupling between intersubband transitions 
in doped quantum wells and LO-phonons. We took into 
account the electron-electron Coulomb interaction and 
we were able to treat exactly the resulting depolarization 
shift. We proved that intersubband polarons can be in 
the ultrastrong coupling regime, reaching extremely high 
values of the coupling constant. We critically discussed 
the relevance of our results both for fundamental and 
applied research. 
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